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Abstract. In this paper we consider a non-smooth vector field Z = (X, Y), 
where X, Y are linear vector fields in dimension 3 and the discontinuity man- 
ifold E of Z is or the usual embedded torus or the unitary sphere at origin. 
We suppose that E is a sliding (stable/unstable) manifold with tangencies, by 
considering X, Y inelastic over E. In each case, we study the tangencies of 
the vector field Z with E and describe the behavior of the trajectories of the 
sliding vector field over E: they are basically closed. 



1. Introduction 

The theory of the non-smooth dynamical systems has been developing very fast 
in the last years, mainly due to its strong relation with branches of applied science, 
such as mechanical and aerospace engineering. Indeed, discontinuous systems are 
in the boundary between mathematics, physics and engineering. 

In many cases, the non-smooth systems are described by systems of ODEs in 
such a way that each system is defined in a region of the phase portrait. The 
boundary of such regions is called the discontinuity manifold [8_ . 

There are a lot of research being made about the local behavior of non-smooth 
systems, near the discontinuity manifold [3]. An interesting phenomenon that oc- 
curs frequently is the sliding motion, when the trajectories on both sides of the 
discontinuity manifold slides over the manifold after the meeting, and there re- 
mains until reaching the boundary of the sliding region j4[ [11] . This phenomenon 
has been studied specially on relay control systems and systems with dry friction 
HIT]. 

The mathematical formalization of the theory of discontinuous systems was made 
precise by Filippov [5] and Teixeira-Sotomayor [121 [TO] , which classified and distin- 
guished three types of regions on the discontinuity manifold (sewing region, sliding 
region and escape region). These regions fully describes what can happen on such 
manifold. Teixeira and Sotomayor |10j have introduced the concept of regulariza- 
tion of discontinuous vector fields, which allows to study phenomena that occurs in 
regular C r vector fields, but in the context of non-smooth vector fields. 

One of the biggest limitations of the theory of non-smooth systems until now 
is the lack of global results |11| . In particular, the discontinuity manifold is often 
taken MS i\ hyperplane. 

This article presents some results for the case where the discontinuity manifold is 
a torus or a sphere in the three-dimensional Euclidean space. We consider inelastic 
vector fields, which are a special class of vector fields with the additional property 
that the discontinuity manifold is composed only of escape and sliding regions 
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(separated by tangency curves) . Inelastic vector fields are widely used in mechanical 
models [2]. 

This paper is outlined as follows: in Section 2, we introduce the concept of non- 
smooth vector fields. In Section 3 we discuss inelastic vector fields. In Section 4 we 
state the main results. In Sections 5 and 6 we proof the main results. 



2. Non-smooth vector fields 

In this section we recall some concepts about non-smooth dynamical systems, 
according to the Filippov convention [S]. We shall adopt M 3 as the ambient space, 
but almost of the definitions and equations can be made for arbitrary dimensions 

Let h : M 3 , — > R, be a function with as regular value and consider the 
manifold M = ^({O}). Also define M+ = /^((O.oo)) and M~ = h^d-oo^)). 
Put A(3) the space of the C-vector fields over K 3 , 0, for r > 2, and 0(3) the space 
of the vector fields Z defined over M 3 as 



Z(x) 



X(x), if h(x) > 0, 
F(x), if h(x) < 0, 



where X, Y £ A(3) and x = (xi, x%, X3). We denote Z = (X,Y). We can identify 
£7(3) = A(3) x A(3) and inherit a topology from this product. 

Now we define Z over M. We break M into three (possibly disconnected) regions, 
according to the following rule (we recall that Xh(x) = X(x) • V/i(x)): 
Mi = {x£ M; Xh(x) ■ Yh(x) > 0} is the sweing region; 
M2 = {x£ M;Xh(x) > 0, Yh(x) < 0} is the escape region; 
M 3 = {x£ M;Xh(x) < 0, Yh(x) > 0} is the sliding region. 

If p G M and Xh(p) = 0, then p is said to be a tangency point for X. If Xh(p) = 
and X 2 h(p) ^ then the tangency is quadratic, and if Xh(p) — X 2 h(p) — but 
X 3 h(p) ^ the tan gency is cubic. 

The trajectories of Z through a point p are defined in the following way: 

a) If p £ M + , then the trajectory of Z through p is trajectory of X through p, until 
the trajectory reach M; 

b) If p £ M~ , then the trajectory of Z through p is trajectory of Y through p, until 
the trajectory reach M; 

c) If p £ Mi, then p connect two regular trajectories of X and Y; we define the 
trajectory through p as the concatenation of these trajectories; 

d) If p £ M2 U M3 then the trajectory of p is the trajectory of the sliding vector 
field associated to X and Y, that is defined as the vector field 

(1) S(X, Y)(x) = ]_ — - (Yh(x)X(x) - Xh(x)Y(x)) . 

We remark that the manifold M is invariant by the flow of S(X, Y). 
c) If p is a tangency point (p £ dM% U dM^) and there are trajectories of the 
sliding vector field starting and ending at p. then the trajectory through p is just 
the concatenation of these trajectories (adapted from Definition 2.5 of [5]). 

So if the closure of the concatenation between a trajectory of the sliding vector 
and a trajectory of the escape vector field is transversal do the curve of tangencies, 
the closure itself is a trajectory. 
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For p in the boundary of the regions Mi, M2, M 3 , the definition of its trajectory 
is much more complicated (or even undefined yet), and is outside our scopus. We 
recommend [6 for these definitions. 

3. Inelastic systems 

In this section we establish some general results on inelastic vector fields. 
Let h : R 3 , — > K., be a function with as regular value and consider the 
manifold M = /i _1 ({0}). Two vector fields X, Y are inelastic over M if 

(2) (Xh)(x) = -(Yh)(x) 

for all x £ M. 

Inelastic vector fields appears naturally in mechanic models |13| and in reversible 
non-smooth models. A planar non-smooth vector field Z = (X, Y) with disconti- 
nuity set E C M 2 is said to be reversible with respect to an involution ip : W 2 — >• M 2 
with Fix(ip) C E, where Fix(ip) = {z 6 I 2 : p(z) = z} if ipX — —Yip. If we take 
ip, E such that E = {(x, y) : x — 0} = Fix(ip) and the vector fields X, Y are even 
in the second variable, then the pair {X, Y) is inelastic over E. 

For x = (x\, X2, X3), we consider ci(x) = x\ + x\ + x 3 — 1 and 02(x) = (x\ + 
x\ + xl + 3) 2 - 16(x 2 + x%). Denote Ej = ^({O}), EJ" = ^((-oo, 0)) and 
E+ = crT ((0, +00)). We note that T,i is a sphere and E 2 is a torus. 

Now we give a characterization of linear E^-inelastic vector fields. 

Lemma 1 (Linear case). Let X,Y be linear vector fields in R 3 . Consider X(x) = 
Ax and Y(yi) = Bx, for A — (a^) and B — (£>ij) matrices. 

a) If X,Y are inelastic over the sphere E 1; then 

/ -ai,i -02,1 - b 2 ,i - 01,2 

B = ^2,1 — a 2,2 

V b 3 ,i h.2 

b) If X,Y are inelastic over the torus E2, then 

(-ai,i -ai, 2 - 02,1 - &2,i -ai,3 
&2,1 — »2,2 — ^2,3 

— &3,1 — a 3,2 —^3,3 

Proof. The equation Q with /i = <7j is equivalent to a linear system of equations. 

Solving this system for the hj's, we obtain the expression of B given in the 
statements. □ 




4. Main results 

For x = (xi,X2,X3), we consider (Ji(x) = x\ + x\ + x\ — 1 and 02 (x) = (x 2 + 
.x 2 + x\ + 3) 2 - 16(x 2 + x\). Denote Ej = ^({O}), EJ = aj 1 ((-00, 0)) and 
Et = ( rTi((0,+oo)). 

Consider X, Y" two linear vector fields defined over an open set !/cK 3 containing 
Ej, j = 1,2. Suppose that X, Y are singular at 0. We can write X(x) = Ax and 
F(x) = Bx, where A = (0^)3x3 and _B = (h hj ) 3x3 . 
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Fix j = 1 or j = 2. Suppose that (X, Y) is a Ej-inelastic pair, that is (Xijj + Yaj)\^ 
0, and consider Zj = (X, Y) the non-smooth vector field 



f X(x) , x 6 E+ 
jl ' ( F(x) , xeET, 



Consider the sliding vector field Sj defined over E* 

(3) gj(x) = j ((y gj )(x)X(x) - (X<Tj)(x)Y(x)) , xeS« 

(Ycx,,-)(x) + (Aer^x) 

where E* = {x € E 3 ; Fa.,- (x) A>,- (x) 7^ 0}. 

Theorem A: If X, Y are linear and inelastic over Ei (sphere), then the sliding 
vector field Si defined on the sphere (if non trivial) has two equilibrium points p± , 
and every other trajectory is closed and equal to the intersection between Ei and 
a plane orthogonal to the line p+p- joining the two equlibria. 

Theorem B: If X, Y are linear and inelastic over E2 (torus) and X<72(x) is qua- 
dratic, then the sliding vector field 5*2 can be defined on S2 except on two circles, 
that consists of curves of singular points. For all other points in E 2 , if non trivial, 
the trajectory of S2 is closed and planar. 



5. Proof of Theorem A 

Before proving Theorem A, we we derive the equations for the tangency points 
and the expression of the sliding vector field. We fix cxi(x) = x\ + x\ + x\ — 1, so 
that Ex = aj~ 1 ({0} is the sphere. 

Let -^(x) = Ax and Y(x) = Bx be two linear vector fields, inelastic over Ei = 
erf 1 ({0}). Note that B is given by Lemma [I] 

Then 

A(7i(x) = 2014X1 + 202^2^2 + 203,3^3 
(4) + (2oi :2 + 2a 2 ,i)a;ia;2 + (2oi :3 + 2a 3i i)xia;3 

+ (2a 2:3 + 2a 3 ,2)x 2 a;3 

According to the classification of quadratic surfaces, due to [Beyer 1987], the 
set {x £ M 3 ; Acti(x) = 0} is or the empty set, or a plane, or the union of two 
intersection planes, or an elliptic cone (the union, by the vertex, of two cones along 
the same axis). So, when not empty, the intersection of {x £ M. 3 ; Atri(x) = 0} 
with the sphere can be or a circle, or two disjoint circles, or two circles intersecting 
transversally. 

Note that if we define Q = A + A 1 , that is, 

(2014 01,2+02,1 ai,3 + a 3j i \ 
01,2 + ^2,1 2o 2 ,2 02,3 + 03,2 
01,3 + 03,1 02,3 + 03,2 2a 3 , 3 J 

then A<7i(x) = x*Qx, so if Q is negative definite, then Xh is the negative of a 
quadratic form, and Acri(x) < for all x. In this case, there are no tangencies. 

Lemma 2. If A is negative definite, then AVi(x) < and X<7i(x) > 0, Vx G S 2 . 
So there are no tangency points over the sphere. 
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Proof. We know by the discussion of the last paragraph that if Q = A + A* is 
positive definite, then there are no tangencies. In |5] is shown that a matrix L is 
positive (negative) definite if, and only if, L + L is positive (negative) definite. 
This complete the proof. □ 

Lemma 3. Suppose that A is a negative definite matrix. Then: 

a) The sphere is a sliding region, and the sliding vector field S2 defined over the 

sphere is given by 

i a 2 ,i + 62,1 03.1 + b 3 ,i \ d 
* 2(X) = I" - 2 V -^^ Z )^ 

«2,1 + &2.1 »3,2 + &3,2 \ 9 



2 2 I 9x2 

'a 3t i + 634 a 3j2 + b 3:2 \ d 



2 2 V J dx 3 

b) The sliding vector field S2 has two symmetric equilibrium points, 

P± = ±-(as,2 + h,2, -(03,1 + &3,l),02,l + &2,l), 
P 

where p = + a\ x + a| 2 + 034 + 634 + b\ 2 + 2b 2 sa 2 ,i + 26 3:2 a 3: 2 + 2034634, 
and all other trajectories of S2 are closed (and contained in the planes orthogonal 
to P-P+ that intersects the sphere). 

Proof. We compute S2 using formula ([3]), and obtain a linear vector field Sa( x ) — 
[Syx, where 





As the matrix [^2] is skew-symmetric, it is clear that zero is an eigenvalue for [S2], so 
there is a line of solution to Sa(x) = 0, that pass through the origin and intercepts 
the sphere exactly in the points p±. The other eigenvalues are a pure imaginary 
number and its complex conjugate. If we consider 

/ 03,2 + 63.2 03,1 + &3,1 , 

( 5 ) v = nr~ ' rx^> 1 

\ a 2,l + "2,1 a 2,l + 02,1 

for a24 + &2,i 7^ 0, then ^(x) J- v, for all x, so the trajectories of ^2 are closed 
(and planar). □ 

Now we back to equation Q and discuss the case where Q (or A, by [5]) is not 
negative definite. 

In this case, following the classification of quadratic surfaces given in [T], we have 
that the solutions over the sphere of the equation X<j\ (x) = describes (see also 
Figure [T]) : 

i) a circle, if the rank of Q is one, 

ii) two disjoint circles, if the rank of Q is two, 

iii) two intersecting (transversal) circles, if the rank of Q is three. 



Then a trajectory of the sliding vector field intersects the tangency set in at 
most 4 points, provided that this trajectory is not a tangency line. In case of finite 
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(a) (b) (c) 

FIGURE 1. Situation where (a) rank(Q) = 1, (b) rank(Q) = 2, (c) 
rank(Q) = 3. The bold lines are tangency circles. 



intersection, recall that from the discussion in the end of Section [2] that we can also 
define the sliding vector field over the entire sphere. 

Now we show that these circles of tangencies are never a trajectory of the sliding 
vector field £2. 

Lemma 4. Let 7 be a non trivial trajectory of S 2 over the sphere. Then 7 intercepts 
the set {x € S 2 ; Xcri(x) = 0} in at most 4 points. 

Proof. The curve 7 is transversal to v (given by Q). There are two options for 
the solutions of Xai(x) = 0, (x, v) = 0: the empty set, or one line passing through 
the origin, or the union of two lines passing through the origin (if 03.2 + 63.2 7^ 0), 
parametrized by 

/,s ( £1. 2&3.1 + £1,203,1 — fl 2,i — &2,l .j. , , \ 
7i,2 (*) = : —r -*,£i,2M > 

V a 3,2 + 3:2 / 

where £i,£2 are the roots of a two degree polynomial S(s) = (given in the Ap- 
pendix). If £1 = £2 then we have just one line of solutions. The lines 7 12 intercepts 
the sphere in none, 2 or 4 points. In particular, the circles of tangencies never are 
a trajectory of the sliding vector field. □ 

So we can consider the (closure of the) sliding vector field on the whole sphere. 
Corollary 5. Lemma^is true even if A is not negative definite. 

6. Proof of Theorem B 
We follow the same sequence of Section [5] 

Fix cr 2 (x) = {x\ + x\ + x\ + 3) 2 - Vo{x{ + so S 2 = ^({O} is the torus. 
Recall that B is given by Lemma [l] 
Then 

Xa 2 (x) = q 2 (x) + g 4 (x) 

where 

92 (x) = —20a lyl x 2 + (— 20a 2 ,i — 20ai j2 ) x 2 xi — 20a 2) 23;| 

+ (—2001,3 + 12a 3j i) X3X1 + (12a 3 , 2 — 20a 2 ,3) x 3 x 2 + 1203,3X3 
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and 

g 4 (x) = 4a M xf + (4a 2 ,i +^ai a )x 2 x\ + (4a M + ^2,2) x 2 2 xj 

+ (4a2,i + 401,2) + 4a 2 ,2a;| + (4ai )3 + 4a3,i) a; 3 a;f 

+ (4a 2 , 3 + 4a3 i2 ) x 3 :r 2 3;i 2 + (4a 1)3 + 4a 34 ) x 3 a;2 2 a;i + (4a 2i3 + 4a 3i2 ) a; 3 x 2 3 

+ (4a M + 4a 3i3 ) x 3 2 xx 2 + (4a 2 ,i + 4a lj2 ) x 3 2 x 2 a; 1 + (4a 2 , 2 + 4a 3:3 ) x 3 2 a; 2 2 

+ (4ai, 3 + 4a 3i i) x 3 3 xi + (4a 2i3 + 4a 3 , 2 ) a; 3 3 a; 2 + 4a 3 , 3 x 3 4 

It is hard to determine the zeroes of Xo~2 for general functions g 2 , q±, So we work 
under the hypothesis 

q 4 (x) = 0. 

In this case, Xa 2 reduces to 

Xcr 2 (x) = q 2 (x) = 32a 3i ia;ix 3 + 32a 3j2 a; 2 x 3 . 

Lemma 6. Under the hypothesis (74(x) = 0, the tangency set over S 2 is the union 
of four circles Ci, C 2) C 3 , C4, where Ci,C 2 are contained in the plane x 3 = and 
C 3 , C4 are contained in the plane 0.3,1X1 + 03 2X2 — (see Figureiity. 



FIGURE 2. Tangency lines over the torus. 



Proof. Just determine the solutions of (72 ( x ) — an d study its intersections with 
the torus S 2 . □ 

Note that we have four intersection points between the circles in the statement 
of Lemma [6] and these curves break the torus E 2 into four regions i?i , i? 2 , i? 3 , R4 , 
considering R\, i? 2 contained in the region x 3 > 0, and i? 3 , i?4 in the region x 3 < 0. 

The expression of the sliding vector field is 

n f s 1/ , \ 9 1, , . d 

62 W = -^(02,1 + 02,1)22 ^ 1- o( a 2,i + °2,i) x i n — , 

2 0x1 2 0x1 

for x G i?i U R 2 U i? 3 U i? 4 - 

The proof of the next two lemmas follow the same techniques of Lemmas [3] and 
[4] and will be omitted. They prove Theorem A. 

Lemma 7. The trajectories of the sliding vector field S2 are transversal to the 
tangency circles C 3 ,C4. So a trajectory of the sliding vector field can pass over 
these circles (crossing regions i?i,i? 2 and R3,Ri). 

Remark 8. Note that if j p is the trajectory of S* 2 passing through p £ S 2 , then 
^2\ ^peRj 7p) = CiUC 2; that is, the curves C±, C 2 are exactly the sets where 

we cannot define S' 2 . 

Lemma 9. The tangency circles C±,C2 are singular tangency points, that is, they 
behave like singular points of S%. 
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Appendix A. Polynomial S(s). 
The polynomial S(s) is given by 

S(s) =Z 2 s 2 + Z 1 s + E , 

where 

S 2 = 202,263,203,2 + 201,163,103,1 + 02,163,263,1 + 02,163,203,1 
+ 02,103, 263,1 + 02,103,2034 + 01,263,263,1 + 01,263,203,1 
+ 01,203,263,1 + 01,203,203,1 + ai,i6 3 ,i 2 + ai,ia 3 ,i 2 + a 2 ,26 3 ,2 2 + a 2 , 2 a 3 ,2 2 

Si = a 3 , 2 3 - 2 62,163,101,1 - 2 62,103,101,1 - 62,163,202,1 - 62,103,202,1 - 62,163,201,2 

- 62,103,201,2 + 63,163,201,3 + 03,163,201,3 + 63,103,201,3 + 03,103,201,3 

- 2 01,163,102,1 - 2 01,103,102,1 - 01,263,202,1 - 01,203,202,1 + 03,163,263,1 
+ 03,103,263,1 + 2 03,263,202,3 — a2,i 2 63,2 — o 2 ,i 2 a 3 , 2 + 03,i 2 6 3 ,2 + o 3 ,i 2 o 3 , 2 

+ 63,2 03,2 + 2 63, 203, 2 2 + 63,2 02,3 + 03,2 2 2 ,3 

So = —01,303,202,1 + 203,363,203,2 + 201,102,162,1 — 02,163,203,1 — 03,163,262,1 

- 02,103,203,1 — 03,103,262,1 — 01,363,202,1 — 01,363,262,1 + Oi,i02,l 2 
+ Oi,i6 2 ,i 2 + 03,363, 2 2 + 03,303, 2 2 - 01,303,262,1 
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